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ABSTRACT OF THE DISSERTATION 


Observation of High Frequency Internal Waves 


in the Pacific Ocean 


by 


Robert Libman Zalkan 


Doctor of Philosophy in Oceanography 
University of California, San Diego, 1968 


Professor Fred N. Spiess, Chairman 


Observations of high frequency internal waves were 
made from 1 to 14 December 1966, in the deep sea off Baja 
California. The fluctuations of the depth of an isotherm 
were measured with a three-element horizontal array 
attached to a stable platforn. 


The waves are characterized as a broad band 
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phenomenon with a continuous distribution in frequency. For 
short intervals of time, a narrow frequency band within the 
continuum is adequately described as a horizontally plane 
wave of a single vertical mode. Furthermore, this simpli- 
fied structure is stationary over time spans of several 
days. High modes are present in the low frequency waves. 
Above 4 cycles per hour, however, the first mode predomi- 
nates, 

The spectral shape is consistent with the shear 
limited equilibrium spectrum proposed by Phillips (1966). 
In addition, the predominance of the first mode at high 
frequencies further emphasizes the importance of shear 
instability in internal wave propagation. 

The horizontal properties of the wave field indi- 
cate well-defined directions of narrow-beamed propagation. 
These directions and the dispersive properties of the 
propagation have led to the identification of local topo- 


graphic features as generating areas of internal waves. 
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CHAPTER I 


INTRODUCTION 


Since internal waves constantly change the interior 
structure and velocity fields of the ocean, they may have 
important effects on dynamic height calculations, current 
measurement, mixing, marine organisms, underwater sound pro- 
pagation, and submarine operations. An adequate description 
of the wave field is required in order to determine the ex- 
tent of these effects. In addition, an understanding of the 
mechanisms of. generation and decay ia necessary to the 
determination of the role internal waves play in energy 
transfer processes. Specifically, the description should 
include modal structure, energy spectra, and directional 
information. 

At the present level of theoretical development, a 
model that adequately describes these features is clearly 
not feasible, and little is known about the mechanisms of 
generation and decay. Our purpose, therefore, is to seek 
such a description from observation of high frequency 
internal waves in the deep sea. 

The complex structure of internal waves makes both 
data collection and analysis difficult. Unlike surface 
waves, internal waves do not have a unique dispersion 


relationship. Instead, there is a different relationship 





for each member of the infinite set of possible modes. In 
addition, the maximum vertical displacements of these waves 
occur at intermediate depths rather than at the free sure- 
face, and are, therefore, difficult to observe. 

To overcome these difficulties, simultaneous measure- 
ments of the fluctuating depth of an isotherm were recorded 
with a three-element horizontal array. This procedure 
allows the establishment of a phase difference (travel time) 
relationship leading to a vector wave number and, therefore, 
directional determination. The wave number magnitudes can 
then be tested against the theoretical dispersion curves for 
confirmation of progressive wave phenomena. 

Stokes (1847) first showed that internal waves can 
exist in a stratified fluid consisting of two homogeneous 
layers. Later Webb (1884) and Greenhill (1887) dealt with 
the multi-layered case. In 19335, Fjeldstad extended the 
theory to include a continuous density gradient. More 
recent systematic treatments of the subject are by Eckart 
(1960) and Tolstoy (1963). A summary of the theory of small 
amplitude internal waves in still water is given in 
Appendix A. 

Several investigators have attempted to obtain di- 
rectional information by towing a thermistor chain in a 
sequence of traverses of different direction. This pro- 
cedure allows one to obtain a directional description of 


the wave field by interpretation of measured Doppler shifts. 





Using this method, Metcalf et al (1962) reported internal 
waves of 17 minute periods traveling south from the equator. 
However, LaFond (1964) in the Pacific and Charnock (1965) 

in the Atlantic, had little success. 

Ufford (1947) appears to be the first to make 
Simultaneous measurements of internal waves at spatially 
separated points. His first method was to make repeated 
lowerings of bathythermographs from the bow, midsection, 
and stern, of a ship located over the San Diego trough, in 
1200 meters of water. The results were inconclusive, due 
to the small phase differences and, therefore, large error 
limits in the determination of wave number. To increase the 
phase differences, Ufford moved to shallow water (40 meters) 
and employed both bathythermograph lowerings from three 
anchored ships and a buoyed array of three thermocouples. 
Unfortunately, each observation spanned only a few hours of 
time, and he used only those portions of the records 
(average 1 hour, 7 minutes) in which a single frequency 
predominated. 

Further work in shallow water (18 meters) has 
been accomplished at the Navy Electronics Laboratory's 
Oceanographic Tower. Here, a three-element array of iso- 
therm followers has been successfully used to determine 
directional information (LaFond 1962a, Cox 1962). Neverthe- 
less, mid-oceanic observations are still required since 


shallow water propagation differs considerably from that in 
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deep water. In pameleaners the wave form,and therefore the 
frequency content,is affected by the nearness of the sea 
floor (LaFond, 1961b) and the propagation direction by shore- 
ward refraction (LaFond,1962b). Also, there is probably a 
preference for the first mode at all frequencies, since the 
shallow depth, and thus narrow waveguide, means large shear 
for modes greater than one. 

We, therefore, seek our information from the deep 
ocean. The difficulty of employing the array technique in 
deep water is attested to by the previous absence of such 


measurements. 





CHAPTER II 
OBSERVATIONAL EQUIPMENT 


A three element horizontal array was obtained by 
attaching heavy duty booms to the Research Platform 
FLIP's (Fisher & Spiess, 1963) bow or upper structure. The 
attachment was done by hinges enabling the booms to be 
properly stowed while FLIP was in the horizontal position. 
Cables to the sea-sensing units were fed over the ends of 
the booms. Upon "flipping" to the vertical position, the 
booms were rigged out to form an equilateral triangle of 
30.5 meters (100 feet) per side (Figure 1). Each boom had 
its own independent sensing system consisting of 1) an 
isotherm follower, 2) a vertical staff of four variable 
frequency thermometers, and 3) a vibrotron pressure guage. 

The isotherm follower is controlled by a thermistor 
in a sensing unit. A servo system is arranged to vary the 
depth of the sensing unit in such a way that the temperature 
measured by the thermistor remains constant. The thermistor 
Signal is developed across'a DC bridge. .One leg of the 
bridge is connected by a cable to a thermistor in the seg- 
sensing element. The other leg is composed of a variable 
bank of resistors for selection of the desired temperature. 
When the thermistor has the same resistance as the variable 


resistors, the bridge is in balance. When the temperature 
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at the thermistor changes, the thermistor changes resistance 
and a DC voltage occurs across the bridge. This voltage is 
fed through magnetic amplifiers and silicon controlled 
rectifiers to the windings of a DC motor. fThe polarity and 
amplitude of the signal determine the direction and speed at 
which the motor turns. The motor drives a winch which 
raises or lowers the sea-sensing element until the bridge 

1s again in balance. In this way the thermistor in the sea- 
sensing unit is kept at the depth of an isotherm and hence 
the name Isotherm Follower (LaFond, 1961la; Good, 1968). 

The upper part of the sea-sensing element is a 
hollow pressure case. A vibrotron pressure guage was instal- 
led in this portion and exposed to sea pressure via an oile 
filled chamber with a flexible diaphragm. The output of the 
vibrotron was used for determining the depth of the sea- 
sensing unit. 

Since servo-systems are subject to dead bands (i.e., 
a certain level of signal is necessary before the system 
responds), a staff with four variable frequency thermometers 
(Squier, 1967) was installed on each sea-sensing unit to be 
used as a vernier. The thermometers were placed .6 meters 
(2 feet) and 1.8 meters (6 feet) above and below the servo- 
system thermistor. A scheme of interpolation was then used 
to determine the actual depth of the isotherm (see Appendix 
B). 


The data were recorded both digitally for later 





analysis and by analog methods for on the spot monitoring. 
The digital data collecting system consisted of 1) a clock, 
2) a sequential sampling unit, 3) a counting system, and 4) 
an incremental tape recorder. Every twenty seconds the 
clock put out a pulse which triggered the sampling unit. 
The sampling unit consisted of a stepping switch which se- 
quentially connected each of the channels of information to 
the counting system for one second'’s duration. The counting 
system was a period counter set to determine the number of 
microseconds required for 4096 cycles of signal. This 
number, along with a serial time count, was then recorded on 
magnetic tape by an incremental recorder. The information 
recorded from each of the sea-sensing units consisted of 
the vibrotron pressure signal, the four variable frequency 
thermometer temperatures, and the vibrotron pressure signal 
again. Each follower was sampled in turn and then the 
azimuthal orientation of the array was recorded to conclude 
the cycle. After the 19 pieces of information had been 
recorded, there was one second of dead time. Then the 
Clock sent out a pulse and the whole procedure started again. 
DC voltages for analog recording were obtained by 
sending the three vibrotron signals and one temperature 
signal to frequency meters. The output voltage from the 
meters drove the pen displacement system of a four channel 
chart recorder, 


since FLIP is free to rotate about the vertical 





axis, it was necessary to employ an azimuthal orientation 
system. An error signal between the desired heading and 
FLIP's actual heading and a rate of rotation signal from the 
syrocompass repeater were used to drive the tilt boxes of 
two hydraulic motors. The motors in turn drove two off-axis 
propellers. The torque from the propellers would rotate 
FLIP until the desired heading was achieved. If FLIP 
drifted off the desired heading the system would again be 
activated to maintain the desired heading. 
A vlock diagram of the data collection and orient- 
ation systems is shown in Figure 2. 
Preliminary data were obtained during February, 
March, and April of 1966 while Flip was operating in 
Hawaiian waters. Analysis of these data indicated the need 
to: 
l. increase the sensitivity and reliability of the iso- 
therm followers 
2. incorporate the thermometer staffs as verniers 
3. change the data recording system to incorporate 2 
4. reposition the booms and winches 
5. change the cable fairleads because of 4 
6. replace lightweight telescoping booms with heavy-duty 
ones 
7. employ the orientation system to keep the array stabil- 
ized in azimuth for as long a time as possible. 


The implementation of these modifications resulted 
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in the system described above. A short seatrial was con- 
ducted near Catalina Island in November. Then FLIP was 
towed to deep water off Baja California and data were 


collected during December 1966. 
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CHAPTER III 
OBSERVATIONS 


A station southwest of San Diego (near 30°N,121°w) 
Was maintained from 1 December 1966 to 14 December 1966. 
FLIP was allowed to drift with the mean integrated current 
over her 300 foot draft. Since the waves travel relative 
to the current, this drift minimized any Doppler effects, 

A track of FLIP's drift is shown in Figure 3. 

Due to the aerodynamic characteristics of FLIP's 
upper structure, she tends to orient with her keel into the 
wind. Once she has self-oriented, the propeller system has 
little trouble maintaining the heading. If the wind 
direction shifts more than twenty degrees, however, the load 
becomes too great on the propeller system and re-orientation 
is required. It was possible, nevertheless, to break the 
data into 15 short series during which the variance of 
array orientation was small. The duration, mean azimuth of 
the keel, and azimuthal standard deviation are given in 
Table I. 

Time records of series 4, 5, and 6 are presented in 
Figure 4, and are typical examples of the internal oscil- 
lations during the entire sampling period. Each series 
contains isotherm depth vs. time from all three followers. 


The middle trace corresponds to the depth scale on the left. 
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Figure 3. Chart of FLIP's track and location 
o of generating areas (after H. W. Menard, 
: “93 1960). Numbers on track refer to the 
day in December 1966 (0800 position). 
Letters refer to generating areas, 
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Division of Data into 15 Series 


TIME (Dec. 1966 


061400-021120 
021120-031005 
021827-031712 
031816-041701 
041701-051547 
051819-061704 
070115-071633 
071633-080625 
0807 51-090636 
090636-100515 
100557-110443 
101430-111315 
120338-122108 
122108-132003 
132003-140556 


TABLE I 


2 STANDARD 
MEAN HEADING (°T) DEVIATION (°T) 


929.2 
239.8 
240.2 
240.7 
240.6 
259.1 
54065 
528.5 
343.2 
342.4 
011.8 
010.4 
002.4 
5S 
007.6 


4.8 
5.2 
5.1 
2.6 
i.e 
evo 
2.7 
ek 
4.3 
1.9 
23 
Diao 
27.35 
6.8 
PORC 
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Figure 4. Depth of the observed isotherm vs time. 


The upper and lower traces have been offset + 10 meters. 
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The upper and lower traces have been offset + 10 meters 
respectively. The interpolation scheme used to calculate 
the isotherm depth may be found in Appendix B. 

The first characteristic one should notice in the 
time records is the predominance of the tidal oscillations. 
Tidal frequencies contain most of the energy and are, there- 
fore, responsible for the largest vertical motions of the 
isothern. 

superimposed on the tidal motions are the high fre- 
quency fluctuations with which this thesis is concerned. 
Several interesting properties of the high frequency 
motions are apparent in the time records. First, the 
records from the three elements exhibit good visual coher- 
ence throughout the data. Second, the larger amplitudes 
of the high frequency oscillations seem to come and go as 
if they were modulated. And finally, the oscillations 
appear to be always present, indicating, at least ina 
qualitative manner, that they may be stationary. These 
conclusions are justified quantitatively in Chapter V. 

The second property, i.e. that of amplitude modu- 
lation, has been observed in shallow water by LaFond (1961) 
and may be caused by two effects. The most likely cause of 
this amplitude fluctuation is the beating effect. When 
energy is present as a continuous spectrum, beating occurs 
among the various frequencies and a rapidly varying ampli- 


tude envelope is developed. 





The second possible cause is a resonance condition 
described by Eckart (1961) and hereafter referred to as 
Eckart Resonance. When the waves are propagated in a well- 
localized group -- localized in space, frequency, and wave- 
number -- a resonant coupling between the layers of high 
stability can occur. The effect of this resonance is to 
cause the energy to oscillate between the layers. If the 
conditions are variable along the propagation path, and one 
measures the motion in one layer only, an amplitude modu- 
lation will be observed, with the modulating frequency 
dependent upon the sharpness of the layers and their verti- 
cal separation. The same effect would be observed if the 
source is not stationary in a spatial sense. That the 
necessary conditions for this resonance existed is shown in 


Chapter IV. 





CHAPTER IV 
ANALYSIS TECHNIQUES 


The power spectrum has become a basic tool for pre- 
senting the properties of long complicated time series in a 
condensed form. Here, the method of first Fourier trans- 
forming and then summing over the squares has been employed. 

Consider a discretely-sampled, real series of finite 
length, X(j), Jj = 0, N=-1. The X(j) may be represented by 


the complex series 


N-| 
A |\277\)m/N dt a a 
X(j)= > X(m) © » J=O,N-I iy 
m=0 
Then the complex Fourier coefficients are: 
A ; N-I \2m7jm/N 
i : = 2 
X(m)=_ 2d X(j) © ,m=0O , N-I (2) 
me 
Taking into account the symmetry of the transform of a real 
series, a spectral estimate of elementary frequency bands 
with two degrees of freedom is obtained by 


S,(m)= Ag, Xm) Rim), mo, M1 a 


where * denotes complex conjunction and A, = 1, m=o 
Am.= 2, m>o 
To increase spectral stability, a sum over adjacent ele- 


mentary bands is formed, 
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Sx (ho) i S,(m) = 5 Am X (m) X(m) (4) 


where p and q are the lower and upper band indices over 
which it is desired to sum and ewes =pt gq. Letting 
&@=q- p+1 (the number of bands included in the sum), 
equation (4) yields 2a degrees of freedom per estimate, 
Since it is desirable to model a continuous spectrum of 
infinite extent which has dimension Lengthé /cycle/time = 
L°t, it is necessary to divide S(h) by the band width, a. 


Thus 


Sx(h) | q A * A 
= = m) X(m 
Eh) = GQ =o 2 Am X (mm) X(m) 5) 
and simple summing is replaced by band averaging. 
The conversion from the indices m, hg, a, to fre- 


quency is accomplished by using the sampling interval At, 
Let NAt = T 


then rin — 2 r(jAt) T 


_N¢ _¢ 
=> ,and Af=7 


so that f= ; uae 


A further discussion of these techniques is pre- 
sented by Tukey in Harris (1966). Haubrich (1965) shows 
that these methods produce the same result as those of 
transforming the correlation function. The actual calcu- 


lations were performed on a CDC 3600 digital computer using 


the fast Fourier transform techniques of Cooley and Tukey 
(1965). 
Once the series has been transformed, the cross 


spectral estimate between two series is simply 


x" 


Gay (h =b5, Aim X (mn) ¥ (mn (6) 


The estimates of coherence and phase difference 


are respectively: 
2 
bee Gy (hg) 
Rye! = E (he) Ey (he) (7) 





and 
yy (Ne) = Arg (Gxy (h,)) y leads x 
| Z (8) 


Considering the triangular array /\ 
3 


» Py(h,) = Arg(G,(h.)), 4=12,23,31 


(9) 


For a single train of plane waves, equation (7) is 
identically equal to 1, and equation (9) is identically 
equal to 0. Loss of coherence and a non-zero Py may be 


caused by a finite beam width (Munk et al, 1963; Cox, 1962). 
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Cox also discusses reduction of coherence caused by turbu- 
lence and variable phase velocity. It should be noted that 
if azl1 (i.e. p = q) equation (7) has only one degree of 
freedom and R (A) is constrained to be equal to one. Sim- 
ilarly, > Pp(h,) in equation (9) is constrained to be equal 
to zero. Therefore it is necessary to band average. 

Failure to average over at least two bands would lead to 

the conclusion of plane wave propagation, when some other 
form of propagation may exist. 

Using the three phase differences defined in 
equation (8), it is possible to determine three vector wave 
numbers. Let vector wave number k have components (ky, Ky). 
The three phase differences, Py for plane wave propagation 


become 


by = Kx Xt + Ky Ye (10) 


where X p and Y) are the spatial separation components be- 
tween elements of the array. Since there are three 
equations and only two unknowns, (kL; Ky) the wave number 
is over determined. Considering the equations two ata 
time, three separate solutions for the vector components 


may be found. 
i/2 


Krms E > | kml] i 


Then 


where 


Kn = ik, “ Jk m 


Ls (12) 
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1 and j being unit vectors in the x and y directions and m 
indexing the three sets of solutions. 


The mean direction of propagation is defined as 


Beam width indication,aéo , is defined as 


AG = |max(Arg(k;))-min(Arg(k;)| i,j=l,3;i#j (4) 
where the branch cut is taken to form Ae as the minimum 


angle between vector solutions k, and k, (as in figure 5). 


Ky 


J 


A@ , 


a k 


Figure 5 
Definition of A g 


It should be noted that if Rg = 1 and, thus,) >, = 0, 
then only two of the equations are independent, the k,'s 
are identical, and“A@=0, These resulta agree with the 


earlier discussion of plane wave propagation. 
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Once a vector wave number with a small A@ has been 
found, it becomes necessary to determine whether the motion 
is due to internal waves or to some other mechanism such as 
anisotropic turbulence. The distinction between wave motion 
and motions due to other causes can be made with the help of 


diagnostic dispersion curves. The solutions,k are, 


rms, 
plotted on an W-K diagram and tested against the theoreti- 
cal dispersion curves to see if there is agreement. It must 
be remembered, however, that solutions which fail the test 
may still be wave motions not due to single plane waves (i.e. 
multiple modes, or two or more plane waves of the same mode 
and frequency traveling in different directions, etc.). 

To obtain the dispersion curves, it is first neces- 
sary to calculate density, P (z), from the observed 
temperature, T(z), and salinity, S(z), distributions. Next 
p(z) is used to calculate N(z) . Finally, for each fre- 
quency band of interest, the wave equation is numerically 
integrated to solve for the k,. 

Figures 6 and 7 show the observed temperature and 
Salinity distributions, the calculated density and stability 
frequency distributions, and the theoretical dispersion 
curves for the first 10 modes. Figure 6 reflects conditions 
just prior to the commencement of wave data collection, and 
Figure 7 reflects the conditions immediately after termi- 
nating data collection. The temperature distribution in the 


upper 200 m was obtained by manually lowering an isotherm 
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Figure 6. Propagation conditions for 1 December 1966. 
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follower. The distributions of salinity and of temperature 
below 200 m were obtained from nearby CalCOFI stations 
(Figure 6 in December, Figure 7 in May). The difference in 
the two salinity distributions at shallow depths is probably 
attributable to sampling at different seasons as opposed to 
being a permanent spatial feature. In any case, temperature 
is the controlling factor within the thermocline and the 


"order of magnitude" 


salinity distributions were used for 
purposes. Since spatial aliasing will occur for > Oo. 
the dispersion curves are terminated at this value. Below 
400 meters T(z) and N(z) monotonically decrease while S(z) 
and © (z) monotonically increase with increasing depth. The 
values of these parameters corresponding to figures 6 and 7 
are tabulated in Table II. 

Several features of Figures 6 and 7 are worth 
further consideration. Both stability frequency distri- 
butions have many relative maxima. These maxima are 
reflected in the dispersion curves as wiggles (Fig. 6, 3rd. 
mode) or changes in curvature (Fig 6, 2nd mode). It is the 
coupling between these layers of high stability which pro- 
duces Eckart Resonance. The greater the vertical 
separation between the layers, the sharper the resonance, 
The resonance regions are shown in the dispersion curves as 
areas where two modal lines converge and then diverge. For 
example, in Figure 7, a resonance region between the second 


and third mode is apparent at a frequency slightly greater 
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than 4cph (\'& 0085). The coupling seems to occur be- 
tween the layers centered about 75 meters and 250 meters. 


Other regions of resonance are also present. 
TABLE II 


Distribution of Temperature, Salinity, Density, 
and Stability Frequency below 400. 
1 December 1966 


Depth Temperature Salinity Density Stability Frequency 


cph) 

500 5.93 34.20 1.02993 2.88 
1000 3.91 34.47 1.03320 2.28 
1500 2.76 34.57 1.03622 2.06 
2000 2.10 34.63 1.03912 1.93 
3000 1.63 34.66 1.04472 Pre 
4000 1.56 34.66 1.05021 1.76 

14 December 1966 

500 5.92 34,22 1.02990 3.00 
1000 3.96 34.49 1.03321 2.28 
1500 2.89 34.58 1.03622 2.06 
2000 2.12 34.64 1.03913 1.93 
3000 165 34,68 1.04474 1.80 
4000 1.58 34.69 1.05024 dealt 


curves show several interesting differences, 


The two sets of stability frequency and dispersion 


The upper 


high stability layer is sharper and slightly deeper in 


Figure 7 relative to figure 6. 


lower layer is weaker and much deeper, increasing the 


On the other hand, the 


separation between layers and causing the 2nd —- 3rd mode 


resonant region to shift to a lower frequency and inverse 


wavenumber. 


While the first modal curve shifted only 
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slightly and only in the low frequency region, all of the 
higher modal curves shifted toward lower frequencies. fhe 
weakening of the intermediate relative maxima may also ex- 
plain why the curves of Figure 7 are much smoother than 
those of Figure 6. Notice, however, that the regions of 
slightly negative curvature still remain in modes 2 and 3. 

These differences may have been caused by both 
temporal and spatial changes in the internal structure. 
But, because the internal waves are continually changing 
the subsurface structure, it is possible that the mean con- 
ditions were stationary throughout the data collection 
interval, and that the differences between figures 6 and 
7 were caused by taking the two profiles when different 
phase relationships existed among the waves. In this case 
these figures would represent perturbations about the mean 
conditions at the two spatially and temporally separated 
points. 

Figure 8 shows the mean depth of the 16.85°C iso- 
therm for series 1-12 and the 16.00°C isotherm for series 
13-15 at each of the three sea-sensing elements. The 
16.85°C isotherm shows large oscillations during the first 
12 series. For the last three series, however, it was 
necessary to shift to a lower temperature in order to re= 
main within the thermocline. Therefore, it is probable that 


the differences are due to a combination of both mechanisms. 
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CHAPTER V 
RESULTS OF ANALYSIS 


The fifteen series tabulated in Table I were 
analyzed by the methods of Chapter IV. As examples, the 
results for series 1, 4, 8, 10, and 14 are shown graphically 
in Figures 10 through 14. The layout of each figure is as 
follows: 

The upper left hand graph shows the three energy 
spectra -— one for each isotherm follower. Spectral averag- 
ing (equation (5) ) was done over eight elementary 
frequency bands yielding 16 degrees of freedom per spectral 
estimate. For 16 degrees of freedom, there is 80% confi- 


2 
dence that 580-<8°=1.47¢ and 95% confidence that Be <s" 


é is the spectral estimate and sé is the 


41.800" where § 
true spectral value. Continuing down the left hand side of 
the page, the second graph shows the three pair coherences 
as defined in equation (7). For 16 degrees of freedom, the 
95% significance level for coherence is .25. The third 
graphs display >) ¢ as defined in equation (9). The upper 
right hand graph shows the computed wave number magnitudes 
(equation (11) ) plotted as points and the theoretical dis- 
persion curves (modes 1 to 10) of 1 December 1966 plotted 


as lines. Since the modal lines converge in the lower fre- 


quency region, an expanded scale plot for series 1 is shown 
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in Figure 9 (corresponds to boxed area of figure 10), # The 
final graph displays the directional information. Each line 
runs from (9-8) to (e+) where © and 46 are defined by 
equations (13) and (14) respectively. 

Various narrow band spectral peaks are apparent 
throughout the 15 series. However, the peak heights lie 
between the 80% and 95% confidence limits and the fre- 
quencies at which they occur change from series to series. 
The peaks are caused by nearby frequencies having the 
proper phase relationships to accentuate their particular 
frequency band during the short duration of sampling. This 
result is consistent with the expected results from sampling 
a continuous, random, Gaussian process. 

The coherence is quite high for the lower frequen- 
cies but in some cases falls off quite rapidly in the 5 to 
7 cph region. Similarly,}? $¢ remains close to zero when the 
coherence is high and then appears to oscillate randomly 
as the coherence drops. The rate at which the coherence 
falls is greatest when the energy spectrum is also falling 
rapidly with frequency and vice versa. A possible expla- 
nation for the fall off is that in these cases there 
appears to be a lack of first mode solutions at the higher 
frequencies. For modes greater than 1, the wavelength 
becomes smaller than the array spacing near these frequen- 
cies. Therefore as frequency increases more and more 


wavelengths occur between the elements of the array (see 
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Cox 1962). In other cases, however, the coherence rises 
again to quite high values and then falls off in the 10 ecph 
region with 5 @$ behaving in a consistent manner. In these 
cases, there appears to be a predominance of first mode 
Bolutions above 5 cph (compare series 1 and 8). 

some of these solutions, particularly above 5 cph, 
lie to the right of the first mode dispersion curve but 
have good coherence, small 2 4, and narrow beam width, Aé®@. 
In several series, these solutions actually seem to form a 
separate modal line below that of the first mode. The best 
example can be found in series 8 where there is a multitude 
of narrow beamed high frequency waves traveling toward 270°. 
The most likely cause of this phenomenon is a Doppler shift 
due to current shear. FLIP drifts with the mean integrated 
current over her entire 300 foot draft and these high fre- 
quency waves are superimposed on the mean integrated current 
over their relatively narrow high stability layer (see Fig- 
ures 6 and 7). The difference between these currents 
required to account for the observed Doppler shifts is on 
the order of only 5 cm/sec. The details of the necessary 
calculations and a specific example may be found in 
Appendix C. 

The solutions considered acceptable as single wave 
propagation, together with their modal numbers and mean 
propagation directions are listed in Table III. Solutions 


falling below mode 1 are included as mode 1 and those 
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TABLE III 


oummary of Acceptable Solutions 


Band No. 2 2 4 5 


PREQUENCY i 9 10 ti i2 eee Gedo. -lreeetcmm—ro 20 21 #22 23 24 25 26 £27 28 29 30 Ba 32 33 34 35 36 on 
CPH mee soo 1. . 
SERIES 05) 1.58 1.74 2.09 2.44 2.79 3.14 3.49 3.84 4920 4.55 4.90 5.25 5.60 5.96 6.31 6.66 7.01 7.36 7.71 8.06 8.42 8.77 9.12 9.47 9.82 10.17 10.52 10.88 11.25 11.50) 110 Ooms 
1 5 1 1 i i 1 i 1 1 1 1 l 1 
Direction 101 082 053 035 o77 102 153118 006 ©55 348 093 O86 112 102 126 108 303 268 332 
Mode 10 4 4 2 R 2 2 1 ir 
Meeectsom 194 207 312 021 165 o47 152 128 of8 111 #2010 083 115 347 136 O75 178 138 150 128 157 
Vode 8 3 3 1 5 ib 2 1 1 i 1 1 il 1 1 1 1 1 1 l 1 il 1 il l 1 
acme 1 1 
2 Direction 340 307 233 039 O75 085 154 148 224 038 078 052 O89 131 138 132 192 161 127 #136 165 231 215 168 254 221 235 220 
, Mode ome 5 3 2 2 2 1 i l i 1 1 1 1 il il ii 1 i i i i i 1 it i 
Direction 233 170 132 036 130 105 010 O45 137 #059 #£=«.023 055 006 020 240 275 286 184 145 249 341 354 185 318 184 280 339 
5 Mode 10 4 5 2 2 5 2 2 1 1 i 1 1 1 1 1 1 1 1 1 1 1 i 1 1 1 
Direction 141 210 095 094 110 114 352 122 069 062 295 121 309 334 315 185 170 166 247 251 276 263 290 225 226 278 
6 Mode 5 3 5 9 5 5 il 2 1 il il 1 1 , 1 1 1 1 l i 1 1 1 1 
Direction 225 046 194 151 233 187 O60 105 154 115 $107 #%121 «330 021 276 263 345 269 257 293 282 oOA6 D7 295 
Mode GeO 1 Oe 10 2 2 2 1 3 2 i, 2 1 l l 1 1 2 2 2 1 i l 1 1 il il 1 i 1 l 
1 Direction 262 322 348 306 194 112 323 232 262 285 019 O12 215 O47 210 315 250 270 240 0352 O59 e522 e299 276 255 261 264 256 279 =+.}276 %315 #42292 
9 Mode non 10 9 © 6 2 3 5 3 2 1 1 2 2 a 1 1 1 i 1 i l l i l il 1 1 i l 1 1 al iL 
Direction 082 301 325 224 027 198 290 344 201 177 4.276 249 261 318 lel eel 233 231 241 313 269 276 268 250 276 271 260 272 277 «2273 293 269 268 272 
Mode 9 10 10 6 1o)6C (UC «4 4 2 i 1 1 2 1 1 1 1 1 il i i a 1 1 1 iL 1 il 1 
Siivirection 317 146 206 276 265 245 275 271 117 #194 (113 214 345 246 225 005 274 265 265 267 255 280 O01 275 241 271 257 271 
Mode lone 10 10 HO: + 29 5 3 2 2 a 2 i a 1 1 il 1 1 i i il a l 
PO Direction 266 260 327 291 358 002 345 051 286 327 Sel 346 336 331 O62 056 227 255 310 282 141 113 010 094 
Mode il 10 10 10 2 5 2 3 2 2 2 iL 1 l 1 1 1 i al 1 il 1 1 1 il 1 il 1 
M A 9 4 5 3 3 2 l 1 2 it i 1 1 1 i i i i i l 1 1 
le ection ae 1 sen 299 277 049 350 270 O11 322 005 278 O06 O54 O24 108 031 065 065 000 147 310 O07 274 339 239 
vy A 3 4 1 2 1 1 1 1 i 1 ql 1 1 1 1 1 i i 1 
a) ton a - ee 349 038 025 358 O41 013 298 036 050 008 017 006 O41 050 OSa ee 002 O07 008 005 
1 1 1 i 2 1 a 1 1 ju es it 1 i 1 1 1 
14 Becton Bes a6 [eo oe ras 015 051 O42 120 315 166 204 319 2e0 190 245 273 285 0235 O42 278 Ov? 
‘ 2 3 6 2 l 1: il 1 2 1 1 i il 1 i 1 1 1 a 1 i 
E> nts on a a aoe ie a ane 047 269 224 O04 322 20 e79 004 035 198 349 245 O87 O34 O31 055 O60 186 353 300 020 
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falling above mode 10 are included as mode 10. All other 
modal indexes refer to the nearest modal line. The 
Solutions were also tested against the 14 December dis- 
persion curves. Complete modal agreement was found except 
at 4.9 cph in series 8, Here, the 1 December curves indi- 
cate third mode and 14 December curves indicate second mode 
but at a end-3rd mode resonant region. 

Although high modes frequently occur at the low 
frequencies, the first mode predominates in the high fre- 
quency range. This feature is probably caused by the wave 
guide narrowing with increasing frequency - thus subjecting 
the higher modes to greater shear until they eventually 
become unstable. In ganeral, it is felt that the modal 
Solutions are quite good. Notice the good agreement between 
the wave number solutions and the theoretical modal lines in 
Figure 9. Credence to this belief is also given by the good 
fit obtained in the slowness plots of the next chapter. It 
must be remembered however that any particular modal So- 
lution, and especially those with large wave numbers, may 
have been shifted relative to the dispersion curves by the 
Doppler effect discussed above. Table IV lists the number 


of Solutions in each mode. 
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TABLE IV 


Number of Acceptable Solutions in Bach Mode 





Mode Number of Solutions Percent of Solutions 
at oT, 66 

2 47 12 

3 23 6 

4 11 3 

2 9 2 

6 5 1 

i 2 1/2 
8 2 1/2 
9 5 1 
10 31 8 
Total D92 100% 


Numerous wave solutions which appear to be station- 
ary over several days may be found by referring to Table III. 
For example the solutions for f, = 2.44 cph show a 2nd mode 
internal wave traveling 153 es 1° during series 1, 2, and 3. 
Similarly, for t. = 7.71, the solutions for series 5, 6, 7, 
8, and 9 show @ = 265 + 15° and for f, = 8.06, 6= 279 = 3° 
during series 8, 9, and 10. 

To define the general shape of the energy spectrun, 


the 45 spectra (three for each series) were ensemble 
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averaged. Considering that the three spectra from 

each series are not independent, and allowing for the over- 
lapping in time of several of the series, the resulting 
spectral estimates have 194 degrees of freedom. For 194 
degrees of freedom there is 80% confidence that 870° <s° 
B30" and 95% confidence that Bic o cer eae A log-log 
plot of the average spectrum is shown in Figure 15. 

The spectrum has been broken into characteristic 
regions according to slope. The slope of the curve and the 
power of frequency at which the spectrum falls are identical. 
Initially (.2 cph to 1.7 cph) where the wave guide (w<N(z)) 
occupies the entire water column, the spectrum falls as 
a, At 1.7 cph, the slope changes to f-1 and the wave 


guide begins to narrow with succeeding spectral estimates 


1.74 eph, Az = 4000 


fee., for tf, for £, = 2.09meculs 


ms; 
Az = 1430 poe and. foReta y= sees acne Az =860,; where Az is 
the thickness of the wave guide). From 5 cph to 12 cph, 
the spectrum falls as pf 77025, Above 12 cph, no internal 
waves appear to be present (the 4 acceptable solutions 
above 12 cph in Table III have been Doppler shifted ), and 
the rate of fall of the spectrum begins to decrease until 
at 50 cph it reaches the inherent noise level of measure- 
ment. This region of zero slope indicates a root mean 
square noise amplitude, C rms atte) of .014 meters. 


The height of the peak at 4.2 cph lies between the 


80% and 95% confidence limits. Also,the individual series 
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solutions for this frequency do not indicate any apparent 
reason for the peak's existence. Therefore, it is felt 
that the peak has no statistical significance. 

The region between 1.7 and 5 cph coincides with the 
region of good wave solutions and the energy appears to be 
attributable to internal waves. Using a Richardson's 
number instability argument, Phillips (1966) has proposed a 
theoretical equilibrium spectrum for first mode propagation 
within the thermocline for frequencies such that the stabil- 
ity layer is small compared to the total depth. His 
spectrum is proportional to f ~1 for kd<< 1 and to f ae 
for kd >> 1 (where d is the thermocline depth). Since pre- 
ceding analysis has shown many solutions of mode higher than 
one within the region kd << ll, the agreement between the ob- 
served spectrum and the equilibrium spectrum indicates that 
the f ~1 power law may apply in general for kd<< 1, and is 
not limited solely to the first mode. The change in slope 
from f ~+ to f7I+25 occurag where there is a rapid narrow- 
ing of the wave guide, and, as Phillips predicted, where 
kd =1 (see Figures 6 and 7). For kd >>1, the first mode 
is predominant and the agreement between the observed 
spectrum ( f 7~2*©5) and the equilibrium spectrum ( f =) 


furnishes additional support to Phillips' Theory. 





CHAPTER VI 
DIRECTIONAL PROPAGATION AND POSSIBLE SOURCES 


To determine whether there were any preferred 
directions of propagation, the accepted solutions (Table 
III) were subjected to further analysis. Each solution was 


spread out over its beam width and normalized so that 
360° 


Sd a (15) 


where S(W)= |/AG : 6@-AG/2 < YS 8+A8/2 
S(p)= O ; w<@-A8/2, p> 8+A89/2 


The solutions were then quantized in one degree intervals 


and added together. 


PU= PSU) 3 j=172, 3/2,...719/2 a6 
| 


where i indexes the solutions and jJ represents the mid- 
direction of the interval. Finally a 10 degree wide run- 


ning mean filter was used to smooth the results. 
+2 


Qi => P(m), f=|,2...360, P(-n)=P(360-n) 7) 


me=)-§ 


The resulting beam patterns are shown in Figure 16, The 
top pattern includes all fifteen series, The remaining 5 
patterns include three series each and show the changes 


which occur with time. Several strong lobes are noticeable 
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Figure 16. Horizontal beam patterns. 


AT 


within the patterns, The prominent lobe toward 270° pro- 
vides an interesting example. It is virtually absent during 
the first 5 series. It begins to appear during the next 3 
series, and for series 7 to 9, overshadows any other 
direction. It then decreases but remains present for the 
final 6 series. The existence of these lobes leads to the 
conclusion that preferential directions indeed do exist. 
The existence of directional propagation leads, in 
turn, to the consideration of possible sources. Azimuthal 
information i8 provided by the beam patterns. To determine 
range, the dispersive character of the wave propagation is 


utilized. The wave energy travels at the group velocity, 


dw _ of 


G, (f)=9 ~ 0X7 (where n indexes the mode) 
which is the slope of the local tangent to the dispersion 
curves at (f,I7An). 22 the source is localized, so that the 
distance of travel is the same for all wave groups, then for 


two wave groups arriving at t and t +At: 


ree Gitar (18) 
X =G(f,) (t+At) (19) 


80 that 


G™ (f,)= 1 a 
XG" (f,)=t+At (21) 


Letting AG = G'(f,)-G"'(f) and subtracting equation (20) 


from equation (21) X= At/AG™ (22) 
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Thus, range to a source is obtained by making slowness plots 
(time of arrival vs. inverse group velocity). The slope and 
teintercept of a reasonable straight line fit give range and 
time of generation respectively. The slowness plots for the 
lobes centered at 270°, 070° and 025° are shown in Figures 
17, 18 and 19. The slope for 270° (2 1.59 day-meters/sec 
gives a distance of approximately 75 miles. Looking in the 
reciprocal direction (i.e., 090°) at a distance of 75 miles 
from FLIP's track, one observes a prominent ridge of sea- 
mounts (labeled "A" in Figure 3). 

In Figure 18 several different slopes are apparent. 
The steepest slope (7.14 day-meters/sec) indicates a distance 
of 3350 miles and again, the source seems to be a cluster of 
seamounts (labeled "B" in Figure 3). The next steepest 
slope (2.18 day-meters/sec) indicates a distance of 100 
miles and the seamount marked "C" in Figure 3 seems to be 
the source. The least slope (.92 day-meter/sec) indicates 
a distance of 43 miles. Since this solution occurs during 
the time when FLIP was rapidly drifting toward the south, it 
is felt that the seamount at "D" (Figure 3) was the source 
and that the wave group also was subjected to advection in a 
southerly direction. 

In Figure 19, only the slope of the heavy line 
(5.00 day-meters/sec) corresponds to a topographic feature 
(range = 233 miles - area labeled "E" in Figure 4). The 


other lines give ranges of 110, 345, and 24 miles. They are 
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Slowness plot for direction 070°. 
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Figure 19. Slowness plot for direction 025°. Numbers refer 
to frequency band numbers of Table III. Super- 
scripts indicate mode. 





Be 


shown to indicate that many of the straight line fits give 
ranges to whith no correlation with topography was found. 
The quantity of solutions used to determine the 
Slopes makesthe range determinations subject to doubt. 
These results were presented not as absolute proof that 
seamounts act as internal wave sources, but merely to 
demonstrate that the solutions point to them as likely pos- 
Sibilities. This conclusion is supported by Sandstrom 
(1966) who has shown that topographic features tend to 


increase the energy in internal waves. 





CHAPTER VII 
CONCLUSIONS 


High frequency internal waves in the deep ocean are 
characterized as a broad band phenomenon with a continuous 
distribution in frequency. The spectral shape is consistent 
with Phillips’ theory of a shear limiting process, and the 
predominance of the first mode at high frequency gives 
additional emphasis to the importance of shear instability 
in internal wave propagation. 

Over short intervals of time, a narrow frequency 
band within the continuum displays great simplicity in 
modal and directional structure and can be adequately de- 
scribed as a horizontally plane wave (i.e., small beam 
width) of a single vertical mode. In addition, it has been 
shown that this simplified structure is statistically 
stationary for time-spans on the order of several days. 

This structural simplicity indicates that the waves are well 
dispersed and implies the existence of fairly localized 
(both in time and space) sources, 

Considering the properties of the total (broad band) 
wave field, additional evidence of localized sources is 
found. The horizontal nature of the wave field displays 
definite anisotropy, indicating well-defined directions of 


narrow-beamed wave propagation. While the exact character 
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of the anisotropy may change appreciably in times of the 
order of several days, well defined directions of narrow- 
beamed propagation always appear to be present. 

The time=~history of the wave field displays the 
presence of dispersive wave trains. As well as indicating 
localized sources, their presence makes possible the calcu- 
lation of source range. These range determinations, coupled 
with the well-defined propagation directions, indicate that 
nearby topographic features such as seamounts act as gener- 
ators of high frequency internal waves, 

In view of the wide distribution of seamounts 
throughout the oceanic basin, the superposition of arrivals 
from the multitude of sources should produce a confused 
"internal sea". Therefore, the waves must not persist over 


long propagation paths. 
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APPENDIX A 


Theory of Small Amplitude Internal Waves 


I. List of Symbols 


C 
f 


a 
x 


Oc CDwTR DPNON <x + as 


[> 
Si 


onan, 
mage” 


velocity of sound 

inertial frequency 

acceleration of gravity 

wave number magnitude 

X=component of wave number vector 
y-component of wave number vector 
mode index 

time 

horizontal coordinate 

horizontal coordinate orthogonal to x 
vertical coordinate (positive up) 
depth of bottom 

force 

stability frequency 

specific volume 

vertical displacement 

density 

angular frequency 

infinitesimal volume 

small vertical displacement 


second derivative with respect to time 


( 9 second derivative with respect to Z 


a 
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II. Infinitesimal Waves 

The following presentation consists of only those 
aspects of internal wave theory necessary to the understand- 
ing of the observational and analytical techniques employed 
in this dissertation. More detailed formulations may be 
found in Eckart (1960) and Tolstoy (1963). 

Consider a vertically stratified, compressible fluid 
at rest ina gravity field. If a parcel of fluid is raised 
vertically from its equilibrium position and then released, 
it will tend to sink with the force: 

2 
F=(p-(p +52 Az) dV)-p 9, AzdV (aa) 
the first term being due to buayancy and the second compres- 
sibility. The vertical equation of motion becomes: 
j= -g2P 2 2dV 
pdaVz= 957 ZdV-p = Z (2A) 
For infinitesimal motions, the solution to equation (2) is 


that of a harmonic oscillator with the frequency: 


N(z)=/-5 55 ~ Poe (3A) 
Equation (34) defines the Brunt-~Vaisala or stability fre- 
quency. In the upper ocean, the compressibility term is 
generally small compared to the buoyancy term and is usually 
neglected. Thus we have harmonic motion caused by buoyancy 


forces due to density stratification. 





ae 


Starting with the three-dimensional equations of 
motion, mass conservation, and invoking the following 
assumptions: 

1) N = N(z) and is independent of horizontal 

coordinates 

2) Mean velocities are equal tc zero 

By), Boussenesq Approximation - Density variation 

may be neglected except where it causes buoyant forces 

4) Nonlinear terms are negligible 


we arrive at the first order equation of motion 


N° (z)-w" = (4 
i we - £2 Vu ie =a) » 


where W is the angular frequency and f is the inertial fre- 
quency due to rotation of the earth and C is the vertical 
displacement from equilibriun. 

If we assume a solution of plane horizontally pro- 


gressive waves of the form 


C(x,y,z,t) =((z) Gee 


(5A) 
and k* =ke ats ky 
then (4A) becomes 
2 2 
, N-w' Wo =0 (6A) 


2 
zz w*-f 
Internal waves have their maximum amplitude below 


the surface of the water column and at this maximum 
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$77 <0 


N*- w? 
ie ee ce wee fe 


C 


Therefore either f <w°<N@ Nese In the ocean N?>f? 


80 that the former condition is the correct one, and f and 
N max are the lowest and highest frequencies permitted to 
free internal waves of form (5A). 
2 2 
If w°>f “then equation (6A) becomes 
2_ 2 
(7A) 


N 
Ct’ ae 


The boundary conditions for our purposes are 
C — O at Jf mm O surface 
C = O at Z2=-D bottom 
The free surface boundary condition may be used at the sur- 
face but has very little effect on the internal modes. 

For arbitrary N(z) and fixed W, solutions to 
equation (7A) are found by iterative numerical techniques - 
starting with C = O at Zz =-D and working to satisfy the sur- 
face boundary condition. For a given W_ , solutions can be 
found for an infinite number of discrete values of K. For 
each such value of fo there is a corresponding amplitude 
function C ne the index n relates to the number of ampli- 
tude maxima within the water column, and (n-1) to the number 
of nodes or zero crossings between the maxima. 


Solutions will be progressive if Gi IN\Z |aomewnene sLaal 
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the water oa (2) will be "trigonometric" where 
n 
w<N(z) and "exponential" where w>N(z) (see Figure 1A). 


N C 





exponential 


trigonometric 


exponential 


ee =: 


Figure 1A. @nd Mode Amplitude Function 


Within the thermocline where the density gradient, 
and therefore N(z), 18 a maximum, the appropriate vari- 
ation of physical parameters and their effects on specific 


volume are given in Table IA 


TABLE IA 
The Variations of Physical Parameters in the 


thermocline and their effects on specific volume 


| da 
Parameter Variation X q 0( ) Specific Volume Change 
Temperature aaa SAxaGaay/ Ac ~1. 6X10"? 
salinity +.3 SE edo Teo ~2,.4x1074 


~4 ~ 4 
Pressure +100db -4X10 /100@b ~4X410 
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It has already been shown that internal waves are oscil- 
lations of isodensity surfaces. Since temperature is the 
dominant parameter within the thermocline, the motions of an 
lsodensity surface may be approximated by those of an iso- 
therm. This approximation is very convenient since 
temperature is much easier to measure than density. Temper- 
ature measurements for this reason have been widely used for 


obtaining observational data. 





APPENDIX B 
Interpolation Scheme 


The following parameters were sequentially sampled 
from each sea-sensing unit: 
Das Ty,» Tos Tz, Ty, Dy 
where D = depth and T = temperature. 


Then the depth D, of each thermometer when sampled 


may be found to first order by 





AD; = (Dp,-D,) 4, i1=+-1,4 
(12) 
Ds = De + AD; +O1 


th 


where on is the vertical offset of the i thermometer from 


the servo-system thermistor. 

The depth of the isotherm (Dp) to first order is 
then found by using the two temperature measurements which 
bracket the isotherm: 

Dp = Dy + | (Dae %% Bid/(T4 =-T42) | (T, = T) 


where T is the temperature of the desired isothern. 
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APPENDIX C 
Doppler Effect 


Phase velocity, c, is defined as 


C=w/k=f/d' 


When sampling at three points simultaneously, the Doppler 
shift occurs only in frequency, leaving the wavelength 


unaffected. Thus: 


Ac = Af/X' 


For frequency 9.47 cph of series 8, \'=8.62xl0°m"! . 


\ '=8.62 x IO °M occurs at 7.90 cph on the first modal 


curve of Figure 6, 


2 
9.47-7.90 y 10° _ 1.57 ~ 5M 


AC= 5 G2x10-3 - 3600. 31 SeC 
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